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In conventional absorption based x-ray computed tomography (CT), the noise variance in recon- 
structed CT images scales with spatial resolution following an inverse cubic relationship. Without 
reconstruction, in x-ray absorption radiography, the noise variance scales as an inverse square with 
spatial resolution. In this letter wc report that while the inverse square relationship holds for dif- 
ferential phase contrast projection imaging, there exists an anomalous scaling law in differential 
phase contrast CT, where the noise variance scales with spatial resolution following an inverse lin- 
ear relationship. The anomalous scaling law is theoretically derived and subsequently validated with 
phantom results from an experimental Talbot-Lau interferometer system. 



DifFcrcntial phase contrast (DPC) imaging and the ex- 
tension to computed tomography (DPC-CT) have at- 
tracted recent interest due to their successful implemen- 
tation using low brilliance sources [1-8]. Initial experi- 
mental results demonstrate that DPC-CT imaging may 
have the potential to quantitatively measure the com- 
position of a material [7] with superior contrast-to-noisc 
ratio (CNR) when compared to conventional absorption 
CT imaging [8, 9]. In order to address whether the sen- 
sitivity of the DPC-CT measurements is sufficient for a 
specific imaging task, two relationships must be studied: 
the noise variance dependence on exposure level and the 
noise variance dependence on spatial resolution. The de- 
termination of first relationship is critical if DPC-CT is 
to be considered for biomedical applications, where the 
minimization of radiation dose is paramount. The sec- 
ond relationship, noise variance vs. spatial resolution, 
will determine whether sufficient CNR can be generated 
at high spatial resolutions and acceptable dose levels. As 
the spatial resolution is increased to improve fine-object 
visualization, the noise variance increases, generally fol- 
lowing an inverse-power law. The power is different for 
projection and tomographic imaging in absorption x-ray 
imaging, and has not been fully investigated in the case 
of DPC and DPC-CT. After determining the above two 
relationships, the potential advantages of DPC and DPC- 
CT over conventional x-ray imaging will become clear, 
allowing for proper selection of applications. 

Recently, the first relationship has been addressed for 
both DPC projection imaging [10] and DPC-CT [8]. It 
was demonstrated that the noise variance in both DPC 
and DPC-CT imaging is inversely proportional to radia- 
tion dose, which is similar to the relationship in conven- 
tional absorption x-ray imaging [11, 12], even though the 
physical mechanism of image formation is dramatically 
different. 

In this letter, we report an anomalous scaling law, 
where the noise variance is inversely proportional to the 
spatial resolution in DPC-CT. This is in stark contrast to 
conventional absorption CT, where noise variance is in- 
versely proportional to the third power of spatial resolu- 



tion [13]. This anomalous behavior indicates that DPC- 
CT imaging may enable higher spatial resolution than 
absorption CT at the same noise variance. In the fol- 
lowing, we will present a theoretical analysis of the re- 
lationship between noise variance and spatial resolution. 
Experimental phantom results are then presented to val- 
idate the theoretical analysis. 

In order to study the noise variance-spatial resolu- 
tion relationship, a Talbot-Lau interferometer DPC and 
DPC-CT data acquisition setup is used. In this setup, 
a partially coherent x-ray beam is diffracted by a phase 
grating with a 7r-phase shift at the mean beam energy. 
There will be a self-image formed at the fractional Talbot 
distance [1, 14, 15], 
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where m = 1,2,3,..., and Zt = 2p^/A is the Talbot 
distance, determined by the wavelength A and the pitch 
p of the phase grating. In order to record the diffracted 
beam intensity modulations, a homodyne technique is 
introduced at the fractional Talbot distance [16]. An 
absorption grating with the same period as the diffracted 
beam modulation pattern {p2 = p/2) is used to analyze 
the beam modulation. After the homodyne analysis, the 
intensity profile is both low-pass filtered and recorded by 
the detector elements. At each detector element {x,y), 
the recorded intensity can be written as [1-4, 10, 17] 
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where; phase shift in the image intensity 0{x. y) is related 
to the phase change $ of the x-ray wave induced by the 
image object by 
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where 5 is the decrement of the refractive index n = 1 — 

5 + il3. Therefore, once phase shifts of the beam intensity 
profiles are measured from different view angles around 
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the image object, an image of the local distribution of the 
refractive index decrement can be reconstructed [3, 18, 
19]. In this letter, ^{x,y) is referred to as the projection 
data and Eq. (3) is referred as the fundamental imaging 
equation of DPC-CT, connecting a measurable quantity, 
(f){x, y), to a line integral. 

In order to measure the projection data 4>{x,y), a 
phase-stepping method [1, 17] is used in which the an- 
alyzer grating is translated by a fraction of the grating 
pitch along the Xg axis: Xg = jp2/M,{j = 1,2, M). 
The measured beam intensity at each phase step is given 

by 
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By multiplying both sides of Eq. (4) by exp(— y27r/M) 
and summing over j, one obtains 
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Namely, after a Fourier transform of the measured in- 
tensity profile over all phase steps, one can obtain the 
desired projection data. 

Due to photon number fluctuations, the projection 
data will fluctuate about a mean value (p, i.e., </> = 0±A(/). 
Using the standard error propagation formula, one can 
calculate that the noise variance ct^ of the projection data 

is determined by the noise variance a'j^j-^ = J^-'^ of the 
measured beam intensity at each phase step: 



where Sij = 1 for i = j and Sij = for i ^ j, and AcpiAc/fj 
is the correlation of the projection data with respect to 
view angle. The indices i,j denote the view angle index 
in the DPC-CT data acquisition. This property is similar 
to conventional absorption CT, where the projection data 
also demonstrates a white noise behavior [11]. 

Using the derived noise model for DPC-CT, we can 
analyze how the noise variance is related to spatial res- 
olution in DPC-CT imaging. DPC-CT images can be 
directly reconstructed using filtered backprojection with 
a Hilbert filtering kernel [4, 18, 19, 21]: 
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where the filtering step F is defined by 
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where lon = 1/(A.t) is the bandwidth of Nyquist fre- 
quency, determined by the target spatial resolution Ax 
of the reconstructed images. In Eq. (11), (p{uj,6) is the 
Fourier transform of the projection data 4>{x, y) at view 
angle 9. Based on this definition, the noise variance of 
the reconstructed image at image pixel {x, y) is 
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Using Eq. (9), one can demonstrate that the Fourier 
transform of the white noise is given by 
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Here the parameter e describes the effective efficiency of 
the interferometer [20]. Prom Eq. (7) one can see that 
the noise variance of the projection data is inversely pro- 
portional to the total mean exposure level MIq at a given 
projection view angle. It is also inversely proportional to 
the square of the efficiency of the interferometer, meaning 
that higher interferometer efficiency will result in lower 
noise variance in the projection data. 

When an image object rotates, the projection data are 
measured from different directions, resulting in data at 
different view angles being uncorrclated. As a result, the 
image noise in the projection data has a white-noise-like 
behavior: 
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where D is the detector element width used for normal- 
ization purposes, and N is the total number of view an- 
gles. Using the above relationship, Eq. (12) can be sim- 
plified to 
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where cr^ is given by Eq. (7). Equation (14) says that the 
noise variance of the DPC-CT image is pixel independent 
and is inversely proportional to the spatial resolution of 
the reconstruction. 

In order to validate the theoretical analysis, we use 
an experimental Talbot-Lau interferometer system con- 
structed at the University of Wisconsin-Madison. The 
data acquisition system consists of three x-ray gratings. 
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(a) (b) 

FIG. 1. Cross-sectional view and reconstruction of the imag- 
ing phantom, (a) Three plastic rods and an air-filled tube 
are inserted into a water-filled cylinder, (b) Reconstructed 
DPC-CT slice from the 2 x 2 binning data set with (80 fim)^ 
voxels. 

a rotating-anode x-ray tube (G1592, Varian Medical Sys- 
tems, California, USA) with a 0.3 mm nominal focal 
spot connected to a generator (Indico 100, CPI, Ontario, 
Canada), a CMOS flat panel x-ray detector (Rad-icon, 
Shad-o-Box 2048, California, USA) with 48 /xm primitive 
detector pitch across a 1024 x 2048 array, and a rotating 
motion stage to enable tomographic acquisitions. The 
three x-ray gratings were fabricated at the University of 
Wisconsin-Madison using similar techniques as described 
in literature [22]. The phase grating was designed to in- 
troduce the differential 7r-phase shift with a 50% duty 
cycle at a mean beam energy of 28 keV. The analyzer 
grating has a pitch p2 of 4.5 /im. In this work, 8 phase 
steps were used, sampled over the 4.5 fj,m period. 

The noise properties were measured in a water-filled 
phantom chamber with an inner diameter of 25.4 mm and 
wall thickness of 1.55 mm. The phantom also contained 
4.76 mm diameter PTFE, PMMA, and POM rods, as 
well as an air-filled tube with an inner diameter of 6.60 
mm and wall thickness of 0.80 mm. A cross-sectional 
view is shown in Fig. 1(a). 

To acquire a complete data set for CT reconstruction, 
360 views of projection data were taken at 1° increments. 
Each projection had a total exposure time of 40 seconds, 
divided over 8 phase steps. The tube potential was 40 
kVp, with a continuous tube current of 20 mA. Once the 
intensity modulation was recorded, the data was pro- 
cessed to extract the differential phase and absorption 
projections. 

Prior to reconstruction, the detector pixels were binned 
1 X 1, 2 X 2, 3 X 3, 4 X 4, and 5 x 5 to explore noise variance 
dependence on spatial resolution. The DPC-CT images 
(Fig. 1(b)) were reconstructed using the FBP algorithm 
outlined in Eq. (10) and (11), while the absorption CT 
images were reconstructed using a standard FBP algo- 
rithm [13]. Two identical scans were performed and the 
reconstructions were subtracted to obtain a noise-only 
image. The subtracted image was divided by to ac- 
count for the additive noise incurred from the subtraction 
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FIG. 2. Log-log plots of relative noise variance against spatial 
resolution, (a) shows the DPC-CT results, where the fitted 
curve is of the form oc Ax'^-^^, with = 0.990. (b) 
shows the absorption CT results, where the fitted curve is of 
the form oc Ax"^■'*^ with = 0.997. Note that the use 
of a log-log plot results in a linear display of the data, where 
the exponent of the fit becomes the slope of the curve. 



of two independent volumes. Due to a geometrical magni- 
fication factor of 1.2, the reconstructed voxel dimensions 
range from (40 /im)^ to (200 /im)"^. 

Figure 2 presents experimental data to demonstrate 
the relative variance dependence on spatial resolution. 
In order to measure the variance, the standard deviation 
was measured in the water background and then divided 
by the mean signal value to determine the relative noise, 
and the relative noise was then squared. The relative 
noise was used to calculate the variance because division 
by the mean signal level makes it a dimensionless quan- 
tity, allowing for comparison across contrast mechanisms. 
The logarithm of the measured data is calculated and fit 
to a linear function against the logarithm of spatial reso- 
lution using a least-squares method. The slope of the re- 
sulting fit is equivalent to the exponential dependence of 
the measured data on spatial resolution. As shown in Fig. 
2(a), for DPC-CT, the noise variance increases with the 
increase of spatial resolution as CTpQQrp oc (Ax)"^"'^, in- 
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dicating that the noise variance is inversely proportional 
to the spatial resolution, as derived in Eq. (14). In con- 
trast, for absorption CT [Fig. 2(b)], the noise variance 
changes with spatial resolution as (7%^^ oc (Aa;)"^'^^. 

In the theoretical derivation, we demonstrated that the 
noise variance in DPC is inversely proportional to photon 
number. Because the number of photons is proportional 
to the area of the detector elements at a fixed photon 
flux, the number of photons scales as the square of the 
detector element dimension. Thus, the noise variance of 
the DPC projection data is proportional to the inverse 
square of detector pitch. It is well known that the same 
property is found in absorption projection imaging [11]. 
Therefore, for projection data, there is no difference be- 
tween absorption and DPC imaging. However, as the 
objective of a tomographic reconstruction is to restore 
the depth information which is lost in projection data 
due to the line integral along the depth direction, one 
may intuitively expect that the noise variance will have 
an additional dimension along the depth direction. This 
naturally results in an inverse cubic dependence on the 
spatial resolution in absorption CT. When the same intu- 
itive argument is applied to DPC-CT, one would expect 
that the same cubic power will appear in DPC-CT, due 
to the fact that both DPC and absorption CT projection 
data have the same noise-spatial resolution dependence. 
Surprisingly, this is not the case, as was shown in the the- 
oretical and experimental results. The noise properties 
of the reconstructions from the two contrast mechanisms 
differ due to the difference in reconstruction algorithms. 
In absorption CT, the ramp filter amplifies the high spa- 
tial frequency content, thus amplifying the noise. How- 
ever, in DPC-CT, the filtering kernel is a Hilbert kernel, 
which equally weights all spatial frequency content. As 
a result, it does not result in an amplification of the high 
spatial frequency noise content. 

Because the main conclusion of the paper is intrinsi- 
cally determined by the filtering kernel of the reconstruc- 
tion algorithm, the same conclusion of the relationship 
between noise and spatial resolution can be drawn for 
diffraction enhanced CT imaging [21] and neutron imag- 
ing [23], as each also measures refraction angle data. The 
data is similarly related to the derivative of a line integral 
as in Eq. (3). 

In conclusion, we have theoretically predicted and ex- 
perimentally validated a novel scaling relationship be- 
tween noise variance and spatial resolution. This scaling 
law dictates that the noise penalty significantly drops at 
high spatial resolution for DPC-CT. For human visual- 
ization, a minimal amount of CNR is often fixed to be 
between three and five [24]. Therefore, because there is a 
smaller noise penalty at high spatial resolution, a smaller 
radiation dose is needed to maintain sufficient CNR for 
visualization. This is an advantage over conventional ab- 



sorption CT, where the dose penalty often hinders the 
application of CT in high spatial resolution imaging. 
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